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I. INTRODUCTION 



Xi-wang Luo, Yong-jian Han, Guang-can Guo, Xingxiang Zho\J3 and Zheng-wei ZhoiJll 

A convenient and effective way in tiie quantum double model to study anyons in a topological 
space with a tensor product structure is to create and braid anyons using ribbon operators connected 
to a common base site [A. Kitaev Ann. Phys. (N.Y.) 303, 2 (2003)]. We show how this scheme can 
be simulated in a physical system by constructing long ribbon operators connected to a base site 
that is placed faraway. We describe how to move and braid anyons using these ribbon operators, 
and how to perform measurement on them. We also give the smallest scale of a system that is 
sufficient for proof-of-principle demonstration of our scheme. 
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' Anyons are exotic quasi-particles in two dimensional systems that obey fractional statistics By the associative 
properties of their underlying algebra, they are divided into Abelian and non-Abelian anyons. Under particle exchange, 
,—1 the wave function of Abelian anyons acquires a phase that can be different than multiples of it. For non-Abelian 
anyons, the wavefunction is subject to a none diagonal unitary gate when two particles are exchanged [H, Anyons 
are important for understanding the physics of two-dimensional strongly correlated systems and they are speculated 
> [ to exist in fractional quantum Hall fluids 
^ ' Aside from their fundamental implications in many-body physics, in recent years anyon systems have been suggested 
^1 as a promising candidate for realizing intrinsically robust quantum computation because of their topological properties 
^H. 0. For instance, Kitaev proposed two exactly solvable spin lattice models with anyonic excitations, the quantum 
I— I' double model and the honeycomb model @, 0|- These schemes are called topological quantum computation (TQC), 

and they take advantage of the topological invariance of anyon qubits to protect against local noise. 
^ ' In spite of the conceptual significance of anyons and their appeal for quantum computation applications, it is very 
I , difficult to study anyons experimentally because the conditions for their existence and observation are extremely 
■ challenging to realize. To date, there has only been some experimental evidence in support of the existence of Abelian 
' anyons in fractional quantum Hall fluids [JQ. For the purpose of universal TQC, Abelian anyons are insufficient 
lO ■ and non-Abelian anyons must be available. There have been several theoretical proposals for direct observation of 
' non-Abelian anyons in quantum Hall systems [1, Q and the honeycomb model [lol - [l^ . Unfortunately, these proposals 
y—( require very complex setup or very low temperatures that are beyond the reach of current experimental capabilities. 

Considering the difficulty in realizing physical systems with genuine anyonic particles, quantum simulation of 
, anyons is very valuable since it provides a viable alternative for studying the kinematics of topological states. In this 
^ ' practice, we will not try to construct the complex many-body Hamiltonian giving rise to genuine anyons which often 
involves interactions between more than two physical particles. Rather, the goal is to find practical methods to create 
topological states in a realizable system and study their properties. Though these states are not truly topologically 
^ ' protected, they can be created, manipulated and observed in a properly designed simulation system, and therefore 
- - are very valuable for the study of anyon physics and topological quantum computing. 

Recognizing the importance of quantum simulation of anyon p hysics, in recent years researchers have proposed 
schemes for simulating both the Abelian [l^ [ij] and non-Abelian (l5l - [T7j quantum double model based on photonic 
and trapped atomic systems. Though these schemes are very valuable for the research of anyon physics, they also have 
serious limitations. In particular, the issue of using topologically protected manipulations only has not been sufficiently 
addressed. For instance, local fluxes and local charges are used in ^16i] to encode quantum information and these local 
degrees of freedom will be disturbed by local noise such as local gauge transformations. In [l^, though information 
in the quantum memory is encoded with anyon states, the computing operations are non-topological. Therefore, 
it is necessary to study new schemes that can demonstrate genuine fractional statistics and simulate topologically 
protected quantum computing. 
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In order to simulate anyon states that can truly demonstrate non-Abelian statistics and universal TQC, we must 
carefully examine the topologically protected space in the system which is the computational Hilbert space for TQC. 
In general, this topologically protected space does not have a tensor product structure and unitary transformations 
of states in it under braiding cannot be described clearly. Consequently, it is unclear how non-Abelian statistics is 
manifested in this space and how it can be used for computation [6,] . To overcome this difficulty, Kitaev suggested that 
one can use an arbitrary but fixed base site as a reference [H to create anyons. By constructing appropriate ribbon 
operators that share the chosen base site as one of their ends (see Fig. [l](a)), one can create anyonic excitations on the 
other end of the ribbons. Since all these anyonic excitations have the same reference base site, an overall topologically 
protected space can be constructed as a tensor product of the topologically protected space associated with each 
ribbon operator @. The corresponding Hilbert space does have a tensor product structure and the transformations 
of anyon states in it under braiding and fusion can be easily derived. Kitaev's scheme offers a feasible and effective 
method for performing topologically protected quantum computation. Unfortunately, a scheme to directly simulate 
this important approach of Kitaev's has not been available, since it is nontrivial to construct appropriate ribbon 
operators required for its realization in a physical system. 

In this work, we show how to dynamically simulate the 6*3 non-Abelian quantum double model by proposing a scheme 
to realize ribbon operators connected to a common base site xq. We describe how to create anyonic excitations using 
these ribbon operators, and how to move the other end of the ribbon to achieve the braiding of anyonic excitations on 
it. We also demonstrate how to detect the braiding and fusion, and propose a method to realize anyonic interferometry 
using controlled ribbon operators to detect the topological states of the anyonic excitations. With these capabilities, 
we can then study the fractional statistics of non-Abelian anyons in the quantum double model, and also simulate 
universal TQC. Based on this, we investigate the minimum scale of a system that is required for proof-of-principle 
demonstration of our scheme, and propose an implementation using superconducting circuits. 



II. THE QUANTUM DOUBLE MODEL 



Kitaev's quantum double model is a spin Hamiltonian which is a sum of quasi-local operators in a two-dimensional 
lattice 0. Its ground states are invariant under gauge transformations generated by some finite group. In this work, 
we consider the non-Abelian group G = S3. For the convenience of discussion and without loss of generality, we focus 
on a square lattice as shown in Fig. [T] Particles live on the edges of the lattice, and their internal Hilbert space is 
described by the group G, V. = {\g) : g £ G}. In Fig. [U arrows are used to label the orientation of the edges in the 
lattice. Reversing the direction of a particular arrow is equivalent to making the basis change \z) — \z~^) (where 
\z) € "H, z € G)for the corresponding qudit. To describe the model, one need to define four types of linear operators 
acting on the Hilbert space H: L^±,g € G and T!^, h E G. Their effect on the basis state \z) is as follows: 

Ll\z)^\gz),L^_\z) = \zg-^), (1) 

T'^\z)^6^,,,\z),T'^\z)^^^-^Jz) . (2) 

To study the effect of these operators on individual qudits, we use j and s to denote an edge of the lattice and 
one of its endpoints. Then one can define an operator Lg{j, s) as follows: if s is the tail of the arrow on edge j then 
Lg(j, s) is acting on the j-th particle. If s is the head of the arrow on edge j, Lg{j, s) is acting on the same 
particle. Similarly, Hp is the left (right) adjacent face of edge j then Th{j,p) is (Tjj!) acting on the j-th particle. 
Local gauge transformations and magnetic charge operators for a vertex s and its adjacent face p are defined as 

Ag{s,p)^Ag{s)= Yl Lg{j,s) (3) 

j^star{s) 

and 

A; 

Bh{s,p)= J2 ]lTh,Ajm,p). (4) 

h\---hk—h m—1 

Here, ji • • • jk are the boundary edges of face p, starting from and ending at vertex s and enumerated in the counter- 
clockwise direction. Notice that Bh{s,p) is a projector into states with magnetic flux h on face p. 
The quantum double model Hamiltonian defined by Kitaev is 



H = Y^{l-A{s))+Y.il-B{p)), 

S p 



(5) 
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where 

e is the identity element of the group. The ground states satisfy A{s)\GS) = \GS}, B{p)\GS) = \GS) for all s and 
p, and exited states involve some violations of these conditions. Because the projection operators A{s) and B(j)) 
are localized, excitations are particle-like living on vertices or faces, or both, where the ground state conditions are 
violated. A combination of a vertex and an adjacent face will be called a site (see Fig. [T|). 

A detailed examination of the excitation properties in the quantum double model is presented in Here, we 
just give a brief description sufficient for our purpose. The quantum double of group G, D{G), is a quasi-triangular 
Hopf algebra that has a set of linear bases Dh^g{x) — Bh{x)Ag{x) on site x — {s,p). Though Dh.g{x) is defined on 
each site, the structure of the quantum double does not depend on the specific site. Quasi-particle excitations in this 
system can be created by ribbon operators F^''^{r) introduced by Kitaev @ which we will discuss in detail later. 
For a system with n quasi-particles, one can use C[n] to denote the quasi-particles' Hilbert space. By investigating 
how local operators Dh,g{x) = Bh{x)Ag{x) act on this Hilbert space, one can define types and subtypes of these 
quasi-particles according to their internal states. It turns out that the types of the quasi-particles have a one-to-one 
correspondence with the irreducible representations of D{G). These representations are labeled by H^|^^ where 

[^] denotes a conjugacy class of G which labels the magnetic charge. R{N^^^) denotes a unitary irrep of the centralizer 
of an arbitrary element in the conjugacy [/i] and it labels the electric charge. The 8 irreps for D{S^) are listed in 
the Appendix. B. Once the types of the quasi-particles are determined they never change. Besides the type, every 
quasi-particle has a local degree of freedom, the subtype. The Hilbert space of n quasi-particles then splits according 
to C[n] — ICdi (S) ■ ■ ■ ® ^rfn ®-^rfi - rfn! where dm is the type of the mth quasi-particle, and ICd,^ is the space 

of local degrees of freedom of the m-th quasi-particle. /Crf,„ is just the representation spaces of (H^'^^ ^s^)m and it is 

spanned by the basis = W(l)) ^ where G [^J]m and \£,(l)) is the basis of the representation 

space of i?(Af[^])m. Mdi- -d„ is the topologically protected space. It is inaccessible by local measurements and 
insensitive to local perturbations. When we braid these quasi-particles, the protected space undergoes some unitary 
transformation, but the type and the subtype of the quasi-particles do not change under braiding. 

Unfortunately, though Mdi - d„ is topologically protected, it dose not have a tensor product structure. This makes 
it difficult to use it directly for quantum computation. However, if we choose a base site and connect it with other sites 
by non-intersecting ribbons (see Fig. [21(a)), we can create quasi-particle excitations whose associated protected space 
can be described with a tensor product structure. Each quasi-particle excitation provides a protected subspace, and 
for a quasi-particle whose type is H^'^^ the corresponding protected subspace is spanned hy {{{ly,^)) = \v) ^ |^)}, 

where i' G [fAm and |^) is the basis of the representation space of R{N[^])m- The overall topologically protected 
space is the tensor product of these protected subspaces. These quasi-particle excitations are non-Abelian anyons. 
Their braiding and fusion rules are given in [Toj , and a brief summary of S3 group and 0(83) anyons is given in 
the Appendix. If the conjugacy class [/i] = [e], where e is the identity group element, the quasi-particle is a pure 
electric excitation, we can simply use |^) — |(j^, 1^)) to denote its topological state. If the unitary irrep R is the 
identity representation, the quasi-particle is a pure magnetic charge excitation, and we use — |(i^, 0) to denote its 
topological state. For these pure charge excitations, their local degrees of freedom can also be simplified in the same 
way. These pure charge excitations based on 0(83) are sufficient for universal TQC [20| . The braiding rules for pure 
charge excitations are 

n\iyi)\iy2) = Wiiy2Vi^)Wi), (7) 

and 

n'W)\C^)^Wm>^)rnnO. (8) 

In the above equations, we have neglected the local degree of freedom. R is the unitary irrep of S3 corresponding to 
pure electric charge excitation, v corresponds to pure magnetic charge, and TZ denotes the counterclockwise exchange 
of the two anyonic excitations. The exchange between pure electric charges is trivial. 



III. SIMULATION OF THE QUANTUM DOUBLE MODEL 



The quantum double Hamiltonian in Eq. ([SJ is very difficult to realize because it requires interactions between more 
than two particles. Instead of engineering a Hamiltonian with genuine multiparticle interactions, we take a dynamic 
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FIG. 1: (Color online)(a) A spin lattice for the quantum double model. The system particles (red) reside on the edges of the lattice. 
All particles on edges that meet at a vertex interact. Ancillary particles (black) reside on the vertices and center of the faces. They 
can interact with their nearest system particles, (b) Short and long ribbons on the lattice indicated by the blue areas. Arrows label the 
orientation of the edges. The short ribbon operator for ribbon ri connects adjacent sites xi and xq and acts on the two red edges. The 
long ribbon connects faraway sites Xi (dashed) and the base site xo- The corresponding ribbon operator acts on the edges in thick lines. 



approach to construct appropriate ribbon operators using local operations involving no more than two qudits at a 
time. 

We focus on the square lattice model shown in Fig. [TJ System qudits live on the edges. In previous schemes [l5l.[l6l| 
based on trapped atoms in 2d optical lattices, ancillary qudits are also used. They are not needed if the system qudits 
are addressable, as is the case in a solid-state system. Nevertheless, in order to demonstrate the wide applicability of 
our scheme, we will present it in a system with ancillary qudits that reside on the vertices and center of faces too (see 
Fig. [T](a)). The ancillary qudits are assumed to be addressable. 

To simulate the quantum double model, we apply single-qudit gates on individual ancillary qudits, as well as single- 
qudit gates on all system qudits simultaneously (since they are not addressable). In addition, we need the two-qudit 
diagonal phase gate U = exjp(i(f)\gi) A{gi \ ® \gi)B{9i\) between an ancillary qudit and its adjacent system qudits. With 
these operations, all two-qudit controlled rotations between the ancilla and system qudits can be constructed efficiently 

MM- 



A. Ground state preparation 



We first need to prepare the system in the ground state, and this can be achieved using a similar method as in 
[l^llB]. Initially, all system qudits and face ancillary qudits are in state |e) and vertex ancillary qudits are in state 
— ^= J2g Iff)- Note that in this initial state, B{p) = 1 for all p. To create the ground state, we only need to make the 

symmetrized gauge transformation projection A{s) Eq. (j6]) on every vertex. This symmetrized gauge transformation 
is carried out from left to right and top to bottom on all vertices in the lattice. The ground state is created as a 
result. Now we take sq in Fig. [l](b) as an example and show how to make the projection A(so). Projections A{s) on 
other vertices are similar. Using the vertex ancillary qudit sq as the controlled qudit, we apply on system qudits the 
controlled gauge transformation J2g \9)so{9\ ® Ag{so). After this operation, the state of system qudit on edge [so, Sq] 
(i.e. the bottom edge of sq in Fig. [T] (b)) is the same as state of ancilla sq. Then, disentangle the ancillary qudits 
from the system by controlled two-qudit gate J2g Iff) [so,s'g] (ff I <8i i+ (sq) between the system qudit on edge [sq, Sq] and 
ancilla Sq. Now the ancilla Sq is in state |e), and we make the operation A{so) — / J2a ^s(*o)- 

Alternatively, one can realize the above disentanglement step by first making proper measurements on the vertex 
ancilla, and then performing appropriate corrections on the system qudits according to the measurement outcome [3l ■ 
For the correction, one needs single qudit unitary operations on system qudits. Assuming one of the operations on 
single system qudit is U , without addressability of system qudits, we can still realize U by preparing an adjacent (and 
addressable) ancillary qudit a in the state \h) and performing the controlled unitary operation |/i)a(^|'Xit/-|-l|G|-i'8'l|G|- 
Similarly, all single-qudit operations on system qudits in our following schemes can be implemented by this technique. 
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B. Anyon creation and braiding 

The central issue in our simulation is anyon creation from ground state. In order to create a topologically protected 
space with a tensor structure that can be used for TQC, we choose an arbitrary site Xq as our base site as shown in 
Fig. [1] (b). The key then is to find ways to create anyons by dynamically performing the ribbon operators F^'^{r) 
connected to ccq on the ground state (see Fig. [1] (b)). This will allow us to create arbitrary topological states of a 
given type by applying the superposition ribbon operator X^^eSa OizF^'^{r)^ where /i corresponds to the anyon type 
and az are the coefficients. 

These ribbon operators F^-^{r) act as follows @ 




(9) 



(yiy2) ^^{yiy2) y's 



The ribbon operators i^(''''f)(r) commute with every projector A{s) and B{p), except when (s,p) is on either end of 
the ribbon. Therefore, the ribbon operator creates excitations on both ends of the ribbon. To simplify our discussion, 
we will use long ribbon operators to create excitations at sites (infinitely) far away from the base site xq. When one 
examines the excitations on these sites, the base site is at an (infinitely) faraway location, and one can ignore the 
effect of the excitations on the base site xq [1]. By doing so, when we apply a long ribbon operator connected to xq 
we can solely focus on the quasi-particle excitation on the other end of the ribbon |6|]. 

In the spirit of this consideration, we will first propose a method to realize short ribbon operators to create anyons 
near the base site Xq, and then show how to move them faraway from xq to implement a long ribbon operator. This 
protocol to move anyons around is needed for anyon braiding and measurement too, and thus is critical ingredient in 
our scheme. 

Now let us study how to perform a short ribbon operator in the general superposition form 

^ a,F^^-(ri), 

where ri is a ribbon connecting the base site xq and a nearby site xi in Fig. [T] (b). This allows us to create an anyon 
at xi in an arbitrary topological state with a certain type determined by /i. 

In order to perform the ribbon operator X^zeSa '^zF^'^ (ri), we need the projection operation '^iO:zi\zi){zi\(zi € 
{z G G\az 0}, i = 1, 2, • • • , to) on system qudit [sq, si] and single qudit gate L'^ on system qudit [si, s[]. Using a 
single projection measurement that corresponds to operation \g){g\, together with appropriate gauge transformations, 
we can realize the projection '^.iCtzi\zi){zi\: 



i 

= azAz^ (^O^sr^ [so,.i] {gMg, \GS) 

i 

= 5^a..A,,(si)|e)[,„,,,](e|A,-i(si)|G5) 

i 
i 

We have used Ag{si) = 1 as there is no excitation at vertex si(i.e. A{si) = 1). Following this idea, we can design the 
following protocol to perform the ribbon operator X^zsSa cizF^'^{ri). 

1. Measure system qudit on edge [sq, si] in the basis {|g)}. This is done by applying the two qudit unitary between 
system qudit on edge [so,si] and ancillary qudit pq, \g)[so.si\{g\ ® L^{po), and then measuring ancilla po 
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in the basis {l^)}- If the outcome is \gi), perform gauge transformation Ag_^{si) (This step is equivalent to a 
projection |e)(e| on edge [so,Si]). 

2. Prepare the ancihary qudit si in the state |0[ct^])sj oc X^i'^^ik') apply the controlled gauge transformation 

Eg\9)sA9\®Msi)- 

3. Measure the ancillary qudit si in the basis {|fc[^]) = Zj^j |0[j,])}, where Z^^^ = '^^ exp{i2TTkj/m)\zj){zj\ and 
\0[z]) — J2i Fo'" outcome fcj^j, apply the single qudit operation Z^^j on the system qudit on edge [sq, si] 
as a correction. Alternatively, this disentanglement step can be done by applying two-qudit gate J^g Iff) [so,si] (ffl® 

(si) between system qudit on edge [so,si] and ancilla si. 

4. Apply operation L^(j, si), where j is the edge [si,s']^]. 

By these steps, we can create an anyon at xi with arbitrary type and topological state. 
We first take the pure magnetic charge excitation 

as an example. Here, is an element of the conjugacy class C. It corresponds to the magnetic flux at site Xi 
and characterizes the local degree of freedom, v, on the other hand, corresponds to the magnetic flux across a 
counterclockwise circle, starting and ending at the base site xq and surrounding only the quasi-particle at xi. It is the 
topological state that we are interested in. Notice that, for this pure magnetic charge anyon i'), the corresponding 

coefficients ^ only when Zi G {z G G\z~^fj,z = v}, and = |C|^ for all Zi £ {z £ G\z~'^fiz — i'}. So, in 
step 2 given above, we prepare the ancillary qudit at si in the state |0[^])sj^ cx jz;), with G {z e G\z~^^z ~ u}. 

As another example, we look at the pure electric charge excitation \£,(l)'-,v) = 1^1^ 

where i? is a irreducible representation of the group G(see the appendix), ^(^^ is the local degree of freedom, and 
T] characterizes the topological state. To create this state, we need to prepare the ancillary qudit at si in the state 
\0[Fi.^ ^ ^])si oc R^^J:^fr|{g^)\gi) lu Step 2. As SL morc generic example, we look at the dyonic combination excitation 

^)) where v and f are both topological degrees of freedom (local degrees of freedom ignored). To create this state, 
we perform corresponding superpositions of ribbon operators X^/ieSa CihF^'^{r) with a fixed g € [v\^ and appropriate 
coefficients at determined by the topological state. 

Once we created anyonic excitations on sites close to the base site xq , we need to move them far away from xo and 
also braid them. For these purposes, we need two basic movements. We need to move the excitation from the original 
site to an adjacent site sharing a common face, and from the original site to an adjacent site sharing a common vertex. 
In Fig. [1] (b), these correspond to moving the excitation from site xi to 12, and from X2 to X3. 

For the first quasi-particle movement from site a;i to site 0:2, mathematically, this means mapping the group element 
corresponding to state for qudit [sq, Si] to a product of group elements corresponding to states for qudits [si, S2] and 
[so,si] (i-e. coherently mapping from \gi)[sa,si]\9'2)[si,s2]\^)rest to |ffi)[so,si]l52)[si,s2]IV'')rest with 5^3^ = gi and |V')rest 
referring to the state of the qudits in the rest of the system) and moving the flux at xi to X2- To do this, we first 
perform the projection operation |e)(e| on the qudit on edge [si, S2](see the first step of anyon creation). Now qudits 
[si,S2] and [sqjSi] are in state |ff2) [^1,^2] l5i)[so,si] with g2= e and g'l = gi, and we have gjffi — gi- The flux at site xi 
doesn't change and now the flux at site X2 is the same as site xi. Then, recall the action of the ribbon operator, we 
need to erase redundant excitation at site xi to finish this movement. This can be done by applying the symmetrized 
gauge transformation A{si) in Eq. at vertex si(for operation A{si), see the ground state preparation). 

The second basic movement from site X2 to X3 has been studied in [l6j . Here we give a similar scheme with a 
simpler disentanglement step. We first coherently map the flux at site X2 to the ancillary qudit at pi by the controlled 
operation 

Mx2,pi) = Bg{x2) (E) L%{pi). 
g&S3 

Then we apply the controlled unitary X^/igSg \^)pi{^\ ® (['S2,S2]) to move the flux from site X2 to site X3, just 
as in |l6']. Finally we disentangle the ancillary qudit pi from the system by first swapping ancilla pi and p2 and 
then applying A{x3,p2)~^ ■ The controlled operation A{x,p) can be decomposed into elementary two-qudit controlled 
rotation operators with each edge surrounding p as a control and the ancilla as the target [l^ . 

The combination of these two basic movements allows us to move any anyonic excitations around. By applying 
them repeatedly, we can realize a long ribbon and move the quasi-particle to an arbitrary site faraway from the base 
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site Xq and also braid the anyons. Notice that, after we create an anyonic excitation at a location faraway from 
the base site Xq, we can use the same procedure to create more anyonic excitations by first applying a short ribbon 
operator and then moving the excitation away. In this process, since the excitations created earlier have been moved 
away from around the base site, they will not affect the short ribbon operations later. 

For TQC based on this quantum double model D{S3), two kinds of ancillary vacuum states are needed, the chargeless 
pair state of pure electric charges |/|_r,|) = } X^j, \v)r ^ \v)r' ^nd the chargeless pair state of pure magnetic charges 

|/[^]) = —/^Y\ ^fEli-i] 1^) ® l'^^^) :20]- Here, 77 and v are topological states and we omit the local degrees of freedom, 

and i? is a irreducible representation of the group S'3(see the appendix). The pair excitations localized at sites Xi and 
Xj(Fig. [l](a)) can be created by a ribbon operator connecting Xi and Xj. If a pair of anyons is created from vacuum, 

the total topological charge of them is zero automatically. So we have |/|_r|) = \R\^ ^Ji{t)C2(i) 

and \I[fj.]) = |C|^ X]z-z-Vi(i)z=A'2(L) "^^^'^'^(^^'■''•'l^'^''' '^'^^'"^ ''i.i ^ ribbon that connects Xi and Xj, and ^i(l), 
^2(L), A*i(L) IJ'2(L) are all local degrees of freedom. Obviously, we can dynamically create this excitation pair 
by performing the corresponding ribbon operators using the method discussed earlier. Therefore, we can simulate 
universal TQC in principle. 



C. Fusion and Topological state Measurements 

We must be able to measure the anyon states in order to simulate the quantum double model. Measurement is also 
needed for TQC. An indirect but technically less challenging method is to perform a fusion measurement. It probes 
the fusion rules of the anyons by which we can confirm the topological nature of the system. If we want to study 
the braiding rules, then a direct measurement of anyons' topological states is necessary. These measurements can be 
accomplished by closed ribbon operator projection and interference. 

There are a few ways to do a fusion measurement. One can first fuse two anyons to one and then measure the type 
of the new anyon. Alternatively, one can move two anyons close and perform appropriate closed ribbon operators that 
encircle them, which does not destroy the two anyons since they are not really fused. In our dynamical simulation, 
this is the realistic approach since we can only move two anyons close rather than fuse them. 

For TQC, the only measurement we need is to detect whether there is a quasi-particle left or whether two anyons 
have vacuum quantum numbers when they fuse [20| . Assuming the closed ribbon c in Fig. [5] (b) encircles only the 



two anyons to be fused, this measurement corresponds to the projection ribbon operator 22 1 



i^E^'-w (10) 

which casts the two anyons into the vacuum quantum number state when they fuse. We then focus on this projection 
ribbon operator and show how to realize it. In principle, projection operators corresponding to other fusion channels 
can be realized in a similar way. 

As shown in Fig. [2] (b), to make the projection in Eg. (1101) . all ancillary qudits are initially prepared in state |e). 
We first prepare ancillary qudit po in state lOrgi) = —i — ^ \g). And if we apply the controlled ribbon operator 

J2geG l5)po(5l ® F3''^{c), we get the final state of the system and ancilla pq: \ipf) oc J2geG \9)po^^''^i^)\^-'i) ■ We can 
then measure po in basis {|fc[G]) = -^[g] l'^[G])}' where Z^^^ = J2gj 6xp(i27rfcj7|G'|)|5j)(gj|. If the outcome is |0[g]), the 
projection in Eq. (|10p succeeds, and the two anyons in ribbon c have vacuum quantum number when they fuse. If the 
outcome is not |0[g]), the projection fails, and there is a quasi-particle left. The controlled closed ribbon operator can 
be realized by first performing a short controlled ribbon operator and then coherently moving one end of the ribbon 
along c until the ribbon is closed. Notice that F^''^(c) contains projections. And for the last two step movements that 
close the ribbon, there may be excitations at the target site of the movements, but our movement scheme above only 
works when there are no excitations at the target sites. So we must carefully close the ribbon. The following steps 
should be followed. 

1. Prepare ancillary qudit po in state |0[gi = —i — ^ \g). Perform the projection operation |e)(e| on the qudit on 
edge [sq, Si] (see the first step of anyon creation). 

2. Apply the controlled two-qudit gate X^gec \9)po{9\ ® Lg(si, [si,s'i\) between the ancillary qudit po and system 
qudit [si, s']^]. 

3. Move the excitation at (si,pi) along the ribbon c in the clockwise direction to site (s„,p„) in Fig. [2] (b). 
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FIG. 2: (a) Non-intersecting ribbons with one common end a;o{the base site). (b)A closed ribbon for fusion measurement. 



4. Measure qudit [sn,so] in basis {\g) : g £ G}. If the outcome is not |e), projection Eq. lfTO)) fails and there is a 
quasi-particle left when the two anyons fuse. If the outcome is |e), perform A{sn) in Eq. ([6]). 

5. Apply the controlled two-qudit gate J2geG l5')po(5'l ® Lg{sQ, [sq, Sq]) between the ancillary qudit po and system 
qudit [sqjSq]. Measure ancillary qudit po in basis |A;[g]) — ^[g]\^[g\)^ where Zj^j = Tig 6xp(*27rA:j/|G'|)|(7j)(gj |. 
If the outcome is IO^qj), the projection is done and the two anyons in ribbon c have vacuum quantum number 
when they fuse. Otherwise, there is a quasi-particle left behind. 

Aside from fusion measurements, we can also measure the topological states of the anyons directly in an interference 
experiment. This can be accomplished by using controlled braiding. In Fig. [T] (b), assume we want to measure the 
anyon at site Xi which is in some unknown topological state \x). We make use of the ancillary anyon at site xi Fig. 
[T] (b) (assuming xi is faraway from the base site xq) which is prepared in the topological state \a). To perform the 
interference measurements, initially, we prepare all ancillary qudits in state |e). Then 

1. prepare the ancillary qudit pi in state + \h)); 

2. coherently braid the anyon at xi around Xi counterclockwise only when the ancillary particle pi is in state \h)] 

3. measure the ancilla pi in basis {"^(|e) i 

The outcome of the measurement on the ancilla pi is m = ±1. The probability distribution for either result is 

P(m = 1) - P(m = -1) = Re[{a\{x\'R^\x)\a)] (11) 

Therefore, changing the state \a) and measuring the probability distributions allows us to obtain the real part of the 
interference amplitudes. Similarly, measuring the ancilla pi in the basis ± *|^))} yields the imaginary part of 

the interference amplitudes. With these interference amplitudes, one can determine the topological state \x) p^ . 

The most critical step for the interference measurement is the controlled braiding. It can be realized by controlled 
movements of anyons and braiding operations. First, we coherently move the anyon at Xi to x^ if and only if the 
ancilla pi is in state | h) . Then we braid anyon at 2:3 around Xi . Finally, we perform a controlled movement of the 
anyon from 2:3 back to xi. 

For the controlled movement from xi to X3 and back, we first perform the projection operation |e)(e| on the system 
qudit [si, S2]. Now the product of group elements corresponding to states of qudits [sq, si] and [si, S2] is the same as 
group element corresponding to state of qudit [sq, si], and the flux at site X2 is the same as site xi. Recall the action 
of ribbon operators, we can make a controlled erase of the quasi-excitations at X2 and xi to finish this controlled 
movement. If the ancilla pi is in state |e), we erase the excitation at site X2 by the symmetrized gauge transformation 
A{s2) Eq. © to cancel the movement. If the ancilla pi is in state \h), we erase the excitation at site Xi by the 
symmetrized gauge transformation A(si) in Eq. ^ to finish the movement. Assume the magnetic flux at site xi is 
which is the local degree of freedom of the ancillary anyon, the controlled movement from site X2 to X3 is done 
by applying the operation |e)pj(e| <S) ^([52,^2]) -I- \h)p-^{h\ (8) L'^{[s2, s'2]) between ancilla pi and system qudit [s2,S2]. 
The controlled movement back is similar. Hence, the controlled braiding can be carried out in the following steps: 

1. Perform the projection operation |e)(e| on the qudit on edge [si,S2]. 
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2. Perform controlled unitary gate |/i)p^(/i| (g) F{si) + \e)p^ (e| (g) between ancilla pi and ancilla si followed by 
controlled unitary gate \e)p^{e\^ F{s2) + \h)p^{h\(^l\G\ between ancilla pi and ancilla S2, where F is the Fourier 
transformation. Apply the controlled gauge transformations J2g lfl')si(ffl ® ^gi^i) ^^'^ J2g 1.9)52(5! ® ^9(^2), 
measure ancillary qudits si and S2 in basis |fc[G])- For outcomes fcs^ and kg^, perform controlled two-qudit gate 

k k 

|e)pj(e| ® Zjgj^ ([si, S2]) + \h)p;^{h\ (g) ([si, S2]) between ancilla pi and system qudit [si,S2] 

3. Perform controlled two-qudit gate |e)pj(e| (E) -^([52,52]) + \h)p-^{h\ (E) L''_^{[s2, S2]) between ancilla pi and system 
qudit [s2, S2]. 

4. Braid the quasi-particle at site 2:3 around site x^. 

5. Apply \e)p^{e\®I{[s2: S2]) + \h)p^{h\i^L''^ {[^2:82]) between ancilla pi and qudit [32,52]: measure qudit [si,S2] in 
basis \g) : g € G. For outcome 51, perform controlled gauge transformation \e)p^ (e|(8)A^-i {s2) + \h)p^ {h\(EAg^ (si), 
then perform ^(§2). 

In order to prevent unwanted braiding caused by noise, we must locate anyons faraway from each other, and rely on 
long braiding operators for interference measurement of the topological states. In contrast, in a fusion measurement 
by closed ribbon projection, when we move two anyons close enough we can choose a very small closed ribbon 
and perform the corresponding projection operator once only. Therefore, generally speaking, the topological state 
interference measurement is much more expensive than the fusion measurement. 



IV. PROOF-OF-PRINCIPLE DEMONSTRATION OF NON-ABELIAN STATISTICS AND TQC 

Since the simulation of non-Abelian anyons is very challenging, we are interested in finding out the smallest system 
that is sufficient for the demonstration of non-Abelian statistics and TQC. Here, we consider the braiding of a pure 
magnetic charge around a pure electric charge, and the fusion of two pure electric charges. And by these two processes, 
braiding and fusion, one could demonstrate non-Abelian statistics and TQC. 

The graph shown in Fig. |3] (a) is the smallest system for demonstration of a unitary TQC gate by anyon braiding. 
We choose site xq as our base site. The gauge transformation Ag(so) at sq is the smallest closed ribbon operator 
around sq. And Ag{so) can also be viewed as creating a pure magnetic charge g, braiding it clockwise around vertex 
So and annihilating it, we denote this braiding process as TZi. If initially there are only pure electric charge anyons in 
the system, then TZi corresponds to a braiding process surrounding the excitation at base site xo (pure electric charge 
anyons live on vertices). Now we consider the process that one create a pure magnetic charge g, counterclockwise 
braid it around all anyonic excitations in the system(excluding the excitation at base site xq) and annihilate it, we 
denote this braiding process as TZ2- So, when initially there are only pure electric charge anyons in the system, 
corresponds to a clockwise braiding process around all excitations in the system(including the excitation at base site 
Xq). And obviously this braiding is trivial, i.e. TZiR-2^ = 1 or equally TZi = Ti-2. 

To demonstrate the non-Abelian anyon braiding, we first prepare the system in the ground state 

|G5)(x \9i')i\9i^)2\gig2'h\g2)A\g2),. (12) 

91 ,92^83 

Then we create a pure electric charge anyon at site a:i, or equally we can say that we create a pure electric charge 
anyon at vertex si (because pure electric charge anyons live on vertices). And we braid a pure magnetic charge anyon 
around this pure electric charge anyon by applying a gauge transformation at sq. Finally we detect the change of the 
anyon state. The detection can be done by the anyonic interference method given above, we make controlled braiding 
of pure magnetic charge anyons around Xi, but we need to repeat the controlled braiding many times to complete the 
detection. Here, in this small-scale system, we can reduce the required number of braiding to just one. Notice that 
the controlled braiding for detection is just the controlled gauge transformation at sq. For example, we create at site 
a;i(i.e. at vertex si) the pure electric charge anyon in topological state 

|0i?2)« Iffr'>l|5r'>2|e)3|ff2)4|52)5 + 

Y Cl.9r^)l|3r^)2|c+)3|g2)4|.g2)5+ (-j^g-) 

E C1,9r'>l|.9r'>2|c_)3|52)4|52)5, 
Si92^^=c- 
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(a) (b) 



FIG. 3: (Color online) (a)The small-scale lattice with a minimum number of 5 system qudits (red[gray in black and white]) for demon- 
stration of the non-Abelian anyon braiding, (b) The small-scale lattice with a minimum number of 7 system qudits (red[gray in black and 
white]) for the demonstration of non-Abelian anyon fusion. 



where i?2 is the two dimensional irreducible representation of ^3 (the group elements and irreps of S3 are given in the 
appendix). Then we apply the gauge transformation At„{sQ), which realizes the braiding unitary gate and transforms 
the anyon at site a;i(i.e. at vertex si) to topological state [lija) due to the braiding rules(see the appendix). To 
detect this change, we prepare the ancillary qudit at sq in state + |c+) + |c_)), apply controlled braiding 

J2g \9)so{9\ ® ^g(so): then measure the ancillary qudit sq in the Fourier basis {\j) = + C"'|c+) + C^"'|c-)),i = 

0, 1, 2} where C = exp(i27r/3). If the outcome is |1), the anyon state is lOfl^). If the outcome is |2), the anyon state is 
l^R^)- Else if the system is in ground state, the outcome must be |0). Note that the states of qudit 1 and qudit 2 are 
always the same in this whole demonstration. If we take the gauge transformation Ag{so) as a single step, then the 
states of qudit 4 and qudit 5 are always the same too. Because of this, we can eliminate one of the duplicates from 
the physical simulator for simplification. We remove qudit 1 and 5. 

For the small scale system in Fig. 131(a), Fig. |l](a) shows the detailed circuit for the described demonstration process. 
We only need one ancillary qudit sq for the measurement, so the total number of qudits for braiding demonstration 
is 4. Initially, all system qudits 2, 3, 4 are prepared in state |e), and ancilla Sq is in state :^(|e) + |c+) -I- |c_)) for 
measurement. For ground state preparation, we need to perform symmetrized gauge transformation A{s) Eq. [6] at sq 
and si. Take A{si) as an example, we can first transform qudit 2 to state lOjSgj) by Fourier transformation, then apply 

controlled unitary |ff)2(ff| €5 i+ (3). To create the anyon at Xi, we first apply the projection |e)(e| on qudit 3. 
Then we need to perform the gauge transformations A^^si) + C^c+ (si) + C*^c_ (si), this is done by first transforming 
qudit 3 to state -^(|e) +C|c+) +C*|c-)) and then performing J2g 15)3(3! 'X' (2). If we do not care about the ground 
state and only want to know the property of the state with anyonic excitations, we could directly prepare the system 
in the excited state as shown in Fig. |4](c) to further simplify our scheme. We can rewrite the topological state [Ofla) 
as 

|0i?.)(x^|.9-^)2|e)3|5)4+ 

Y,C\g-'c^%\c+h\g)4+ 

a 

J2C\9-'cZ')2\c-h\gU. 
a 

To create this state, we first transform the system to state 

i=|e>2(|e>3 + C|c+>3 + C1c_)3)|0[53]>4 

by single qudit unitary, where \0[S3]} — |e) + |c+) + |c_) + \to) + \ti) + |t2). Then apply the controlled unitary 
Eg \9)4{g\ ^ L<L{2) followed by Eg 15)3(5! ® ^-(2). This is shown in Fig. H (c). 

For the fusion measurement demonstration, we need more sites to place excitations. The graph shown in Fig. [3] 
(b) is the smallest system for demonstration of this elementary operation, where Xq is the base site. We only consider 
the fusion measurement of pure electric charge excitations. We first prepare the system to ground state, then create 
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FIG. 4: The detailed quantum circuits for the demonstration of anyon brading and measurement. The Fourier transformation F 
transforms state |e) to lOfg^j). The controlled Ug = Ylg \9) control{g\ ^ Ug(target). Ui is an unitary operation that maps state |e) to 
"^(k) + C|c+) + C*|c-)),C = exp(i27r/3). M, Mi and M2 are measurements in basis {|fe[z]) = •^[z]|0[z])>-2^ = SjC+jC-} and 

= -Zj^^j |0[53] )}. (a) is the quantum circuit for ground state preparation, anyon creation and braiding, and detection, (b) is the 
quantum circuit for ground state preparation, anyon creation and fusion detection, (c) and (d) are the corresponding quantum circuits 
with direct anyon state preparation. 
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FIG. 5: (Color online)(a) A photonic qudit based on six TLRs that are capacitively coupled to current biased Josephson junction. (b) 
A 4-junction SQUID device is used to interact photonic qudits. (c)The superconducting circuit for demonstration of anyons braiding 
statistics, (d) The superconducting circuit for demonstration of the fusion of anyons. System qudits {red[gray in black and white] square) 
and ancillary qudits (black square) are coupled by the 4-junction SQUID (blue[gray in black and white] circle). 



two arbitrary pure electric charge anyons at sites xi and X2(i.e. at vertices si and §2). Since there are only the two 
anyons in the system, and the total topological charge of all excitations in the systeni(including the excitation at base 
site) is vacuum. So iff there is no excitation at the base site xq, or equivalently iff the state of the system lip) satisfies 
A{so)\ip) = \ip){al\ excitations have no magnetic charge, so B{p) = 1 for all p), the two pure electric charge anyons at 
si and S2 fuse to vacuum. So to detect the fusion, we prepare the ancillary qudit at sq in state |0[53]), apply controlled 
gauge transformation l5)so(ffl ^gi^o), then measure the ancillary qudit sq in the basis |fe[S3]) = Zj^^j jOjSg]). If 
the outcome is "lOfg^])", we get vacuum after fusion. Otherwise, we get a quasi-excitation left behind after fusion. 

Since the states of qudit 1 and 2 are always the same in the whole process, so are those of qudit 6 and 7, we can 
eliminate qudits 1 and 7. Also we need an ancillary qudit sq for measurement, so the total number of qudits for fusion 
demonstration is 6. As an example, we can create two electric charge anyons in topological state \0r2)si\0r2)s2^ and 
measure their fusion. Based on the same consideration as the braiding demonstration, shown in Fig. 2] (b) and Fig. 2] 
(d) are the quantum circuits for this process with and without the ground state preparation respectively. Notice that 
the two anyons at si and S2 are in the same state 10^,,), the operations on qudits 2,3 are the same as operations on 
6,5 due to this symmetry. 

Beside single qudit unitary operations on individual qudits, the above demonstration processes need two qudit gates 
between adjacent system qudits i and i + For measurement, we need two qudit gates between ancilla sq and system 
qudits 3, 4, 5. 

For physical implementation, we can simulate this model using a variety of physical systems such as cold atoms or 
trapped ions. Here, we point out the interesting possibility of simulating it on a fully integrated superconducting chip 
which provides a powerful platform for manipulating and interacting photonic qubits based on on-chip transmission 
line (TLR) resonators [2^. We briefly describe the system only and leave the details to elsewhere. As shown in 
Fig. [S] [23], 6 transmission line resonators are used as 1 qudit. There is only 1 photon in the system of these 6 
resonators and its location represents the 6 internal states of the qudit. By coupling the TLRs capacitively to current 
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biased Josephson junctions (CBJJ) as shown in Fig. [5] (a), we can shift the relative energies of the TLR modes and 
exchange photons between adjacent TLRs |23)], and thus reahze arbitrary single qudit gates. Further, we can use 
a four junction SQUID shown in Fi g, p (b) to interact photons in 2 adjacent qudits to realize the 2-qudit phase 
gate U = exp{i(j)\gi)A{gi\ <8) \gi)B{gi\) [23j . Thus, we can perforin all operations required to simulate the non-Abelian 
anyons in this system. Shown in Fig. [5] (c) and (d) are possible designs of the smallest superconducting circuits for 
the demonstration of anyon braiding and fusion measurement. 

V. CONCLUSION 

In summary, based on Kitaev's quantum double spin lattice model, we have proposed a method to simulate the 
non-Abelian statistics and universal TQC by performing appropriate ribbon operators. In contrast to earlier studies, 
our scheme is based on braiding of anyons created by long ribbon operators connected to a common base site, and 
hence can simulate genuine anyon states in a topologically protected global space. By designing the smallest system 
sufhcient for the demonstration of fractional statistics and TQC, we have shown that the requirement of our proposal 
is only modest and it is an attractive scheme for experimental studies. 
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Appendix. A. 

Some properties of the groupS'3 : 

5*3, the group of permutations of three objects that we label 0,1,2, is the smallest non-Abelian group. S3 contains 
three conjugacy classes, namely: 

1. Identity e. Centralizer A^^gj = S3. 

2. Reflections to = (01), h = (12), t2 = (20). Centralizer N[t] = {e,to} = ^2. 

3. 3-rotations c+ — (012), c_ = (021). Centralizer N^^] = {e, c+,c_} = Z3. 
The multiplication rules for S3 are as follows: 

tiU = e,tjtk = c^^^Jorj ^ k. 

tiC± = ti±l,C±ti — tizpl, 

CpCp ^ c^p, CpC„ = e, J or a ^ p. 

The operations involving e are trivial. Here ej^k — ± for k — j ± 1 (modulo 3) respectively. 

The group has three irreducible representations (irreps). Two one-dimensional irreps are the trivial one 
and the signature representation Ri(e) — Ri{cp) = -f 1, Ri{ti) = — 1. The two dimensional representations are 

2n 27r 
i?2(e) = l2,R2{tk) = a''exp{i — ka^),R2{c±) = exp{±i—a^). 

Explicitly, 

i?2(e) = ( J ? ) '^2(c+) = ( e ) '^^^'-^ = ( C ) ' 
-^2(^0) =(^10) '^2(^1) ^ ^ ^0 ) '-^2(^2) ^ {^^* ) ' 
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where ^ = e*^'^/^. 

The characters Xr{9) = i^[R{9)\ ^re equal to ±1 for one dimensional reps. For two dimensional reps, 

Xfl2(e) = 2,XR2{U) = 0,Xi?.2(cp) = -1- 

The permutation representation for 5*3 is a set of 6 x 6 matrices that faithfully represents group left action on the 
basis {|e), \to), |ti), jtz), |c+), |c_)}, i.e. L'l\g) = \hg). Similarly the right multiplication L'i|g) = Igh''^), [L'1,L'^] = 
and the unitary matrices are given by 
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Further, the group has a semi-direct product structure which may be exploited to simplify physical realizations. 
We have 5*3 = Z3 Z2 = {a, h\a^ = e, 6^ = e, : cj)b{a) = bab^^ — 0?}. We can choose Z3 = {e, c+, c_}, Z2 = {e, to}, 
and any element g £ S3 can be written as g = c^tQ for r € {0, 1, 2}, s € {0, 1}. The above unitary matrices can be 
simplified in {|7')|s) = |c!j_ig)}, a product basis for a qutrit and qubit [Tg} . 



15 



^+ 


— -L3 i^y -L2 J 


— F(^ 91 6?) rr^ 


^ -\- — ^ ^ } 5 


^ + 


— r \y, Ij 'i9 U 




-1-2, ^-f- — (X) i2j 


LI 


-l3(»l2,i*° 


= l3«)Cr"=, 




L'l 






: X^^ (g) (T+ + X (g) <T 




=x-i® |0)(0 








=X® |0)(0| 4 


|1)(1|, 





where F{i,j) = (|«)0'| + ffl 1 flip two basis states of qutrit, and X = J2i N + 

Appendix. B. 

Some properties of the 0(83) anyons : 

The types of the quasi-particle of the Quantum Double model have a 1-to-l correspondence with irreducible repre- 
sentations of D{G). The 8 irreps for D{S3) and the corresponding quantum dimensions are: 

• Vacuum IlK , d = 1; 

^1 

• Pure electric charges n'*^' , nfe' ,d~l,2; 

-R;^ ' -'1^2 ' ' ' 

• Pure magnetic charges n|j^ , n!y*] , d = 2, 3; 

• Dyonic combinations nj^j , n^_| , Illy*] . c? = 2, 2, 3, 

where 70,71 correspond to the identity and signature representation of iV^^] = {e,io} — ^2- PoiPi,P2 are the three 
one dimensional irreps of N\^c] — {e, c+, c_} = Z3. /3o is the identity representation. Here we only give braiding rules 
of the pure charge excitations, a full description can be found in [6l. Il9l|. 

where R is the unitary irrep of S3 and corresponding to pure electric charge excitation, ly corresponds to pure 
magnetic charge, TZ represents the counterclockwise exchange of the two anyonic excitations. The exchange between 
pure electric charges is trivial. 
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